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Gaussian matrix product states are obtained as the outputs of projection operations from an ancillary space 
of M infinitely entangled bonds connecting neighboring sites, applied at each of TV sites of an harmonic chain. 
Replacing the projections by associated Gaussian states, the building blocks, we show that the entanglement 
range in translationally-invariant Gaussian matrix product states depends on how entangled the building blocks 
are. In particular, infinite entanglement in the building blocks produces fully symmetric Gaussian states with 
maximum entanglement range. From their peculiar properties of entanglement sharing, a basic difference with 
spin chains is revealed: Gaussian matrix product states can possess unlimited, long-range entanglement even 
with minimum number of ancillary bonds (M=l). Finally we discuss how these states can be experimentally 
engineered from A'^ copies of a three-mode building block and TV two-mode finitely squeezed states. 
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Introduction. — The description of many-body systems and 
the understanding of multiparticle entanglement are among 
the hardest challenges of quantum physics. The two issues 
are entwined: recently, the basic tools of quantum informa- 
tion theory have found useful applications in condensed mat- 
ter physics. In particular, the formalism of matrix product 
states (MPS) fT] has led to an efficient simulation of many- 
body spin Hamiltonians |2] and to a deeper understanding of 
quantum phase transitions |3]. 

Beyond qubits, and discrete-variable systems in general, a 
growing interest is being witnessed in the theoretical and ex- 
perimental applications of so-called continuous variable (CV) 
systems, such as ultracold atoms or light modes, to quan- 
tum information and communication processings |4]. Besides 
their usefulness in feasible implementations, quasi-free states 
of harmonic lattices, best known as Gaussian states, are en- 
dowed with structural properties that make the characteriza- 
tion of their entanglement amenable to an analytical analysis 
fs']. Within this context, the extension of the matrix product 
framework to Gaussian states of CV systems has been recently 
introduced as a possible tool to prove an entropic area law for 
critical bosonic systems on harmonic lattices |6], thus com- 
plementing the known results for the non-critical case |7]. 

In this work we adopt a novel point of view, aimed to com- 
prehend the correlation picture of the considered many-body 
systems from the physical structure of the underlying MPS 
framework. In the case of harmonic lattices, we demonstrate 
that the quantum correlation length (the maximum distance 
between pairwise entangled sites) of translationally invariant 
Gaussian MPS is determined by the amount of entanglement 
encoded in a smaller structure, the 'building block', which 
is a Gaussian state isomorphic to the MPS projector at each 
site. This connection provides a series of necessary and suf- 
ficient conditions for bipartite entanglement of distant pair of 
modes in Gaussian MPS depending on the parameters of the 
building block, as explicitly shown for a six-mode harmonic 
ring. For any size of the ring we show remarkably that, when 



single ancillary bonds connect neighboring sites, an infinite 
entanglement in the building block leads to fully symmetric 
(permutation-invariant) Gaussian MPS where each individual 
mode is equally entangled with any other, independently of 
the distance. As the block entropy of these states can diverge 
for any bipartition of the ring |9], our results unveal a basic 
difference with finite-dimensional MPS, whose entanglement 
is limited by the bond dimensionality |10| and is typically 
short-ranged 111). Finally, we demonstrate how to experimen- 
tally implement the MPS construction to produce multimode 
Gaussian states useful for CV communication networks 1 12]. 

Gaussian matrix product states. — In a CV system consist- 
ing of TV canonical bosonic modes, described by the vector 
R = {^ij 92, • • ■ , qN,Pi,P2 ■ ■ ■ ,Pn} of the field quadrature 
operators, Gaussian states (such as coherent and squeezed 
states) are fully characterized by the first statistical moments 
(arbitrarily adjustable by local unitaries: we will set them to 
zero) and by the 2TV x 2TV real symmetric covariance matrix 
(CM) 7 of the second moments 7^ = l/2{{Ri, Rj}) 01. 

The Gaussian matrix product states introduced in Ref. 0] 
are TV-mode states obtained by taking a fixed number, M, of 
infinitely entangled ancillary bonds (EPR pairs) shared by ad- 
jacent sites, and applying an arbitrary 2M 1 Gaussian op- 
eration on each site i = 1, . . . , TV. Here the cardinality M 
is the CV counterpart of the dimension D of the matrices in 
standard MPS ||ljj] . Such a process can be better understood by 
resorting to the Jamiolkowski isomorphism between (Gaus- 
sian) operations and (Gaussian) states |8]. In this framework, 
one starts with a chain of TV Gaussian states of 2M+1 modes: 
the building blocks. The global Gaussian state of the chain is 
described by a CM T = 0,^1 7^*'- As the interest in MPS 
lies mainly in their connections with ground states of Hamil- 
tonians invariant under translation |6], we can focus on pure 
(Dct7['l = 1), translationally invariant (7['l = 7 Vi) Gaus- 
sian MPS. Moreover, in this work we consider single-bonded 
MPS, i.e. with T\f = 1. However, our analysis easily general- 
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FIG. 1: (Color online) Gaussian matrix product states. F'" is the 
state of A'^ EPR bonds and 7 is the three-mode building block. After 
the EPR measurements (depicted as curly brackets), the modes 7^ 
collapse into a Gaussian MPS with global state r°"*. 



izes to multiple bonds, and to mixed Gaussian states as well. 

Under the considered prescriptions, the building block 7 is 
a pure Gaussian state of the three modes with respective CM 
£*i,2,3- As we aim to construct a translationally invariant state, 
it is convenient to consider a 7 whose first two modes have 
the same reduced CM. This yields a bisymmetric |9], pure, 
three-mode Gaussian building block whose CM 7 can be writ- 
ten without loss of generality in standard form lITsll . with 
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s and (Det 0:3)^/^ = x. This 



choice of the building block is physically motivated by the fact 
that, among all pure three-mode Gaussian states, bisymmetric 
states maximize the genuine tripartite entanglement 113(1 . It is 
instructive to write 7 in the block form 
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where 7^^ is the CM of modes 1 and 2, 7^ is the CM 
of mode 3, and the intermodal correlations are encoded in 
Isx- Explicitly ill: 7s^ ^ *+ ^ 
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The MPS construction works as follows (see Fig.[T). The 
global CM r ~ 7 corresponds to the projector from the 

state r'" of the N ancillary EPR pairs, to the final A^-mode 
Gaussian MPS r°"*. This is realized by collapsing the state 
r*", transposed in phase space, with the input port F^s = 
- 7j^ of r, so that the output port Fj; ~ - 7^ turns into 
the desired F°"*. Here collapsing means that, at each site, 
the two two-mode states, each constituted by one mode (1 or 
2) of 7j,^ and one half of the EPR bond between site i and 
its neighbor (i — 1 or i + 1, respectively), undergo an "EPR 
measurement" i.e. are projected onto the infinitely entangled 
EPR state BH]. An EPR pair is described by a CM 
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which corresponds to a two-mode squeezed state of modes i 
and j, in the limit of infinite squeezing (r — * cxo). The input 



state is then F'" = limr^oo 0f'^ ""i/i+i (j'), where we have 
set periodic boundary conditions so that + 1 = 1 in labeling 
the sites. The projection corresponds mathematically to taking 
a Schur complement (see Refs. |6, 8] for details), yielding an 
output pure Gaussian MPS of N modes on a ring with a CM 



pout pT (p^^ ^ 0r"0)-iF, 



(3) 



where F^x 0^ 7^^, and = 0^ diag{l, 1, -1, -1} rep- 
resents transposition in phase space {qi qi, pi ^ —Pi)- 

First we note that the Gaussian states constructed in this 
way are ground states of harmonic Hamiltonians (a property 
of all Gaussian MPS |6]). This follows as no mutual corre- 
lations are created between the operators qi and Pj for any 
i, j = 1, . . . , A^, having both EPR bonds and building blocks 
in standard form. The final CM Eq. Q thus takes the form 



po«t ^ (7-1 . 



(4) 



where C is a circulant N x N matrix. It can be shown that a 
CM of the form Eq. corresponds to the ground state of the 
quadratic Hamiltonian H = ^{J^iPl + Y^ij QiVij<lj)^ with 
the potential matrix given hy V = fl4ll . 
Entanglement distribution. — In the Jamiolkowski picture 
j^l^, different MPS projectors correspond to differently en- 
tangled Gaussian building blocks. Let us recall that, accord- 
ing to the "positivity of partial transposition" (PPT) criterion, 
a Gaussian state is separable (with respect to a 1 x A^ bi- 
partition) if and only if the partially transposed CM satis- 
fies the uncertainty principle |15]. As a measure of entan- 
glement, for two-mode symmetric Gaussian states 7^ ^ the 
entanglement of formation Ep is computable via the for- 
mula |16]: Epifij) = max{0, /(7/ij)}, with /(x) 

(I+X)^ , (l + xf (l-xf 1 (1-x)^ XT -t- 

4^ log - log ■ Here the positive pa- 

rameter rjij is the smallest symplectic eigenvalue of the partial 
transpose of 7, ^ . For a two-mode state, rjij can be computed 
from the symplectic invariants of the state flTil . and the PPT 
criterion simply yields 7, ^ entangled as soon as rji^j < 1, 
while infinite entanglement is reached for rji_j — > 0+. 

We are interested in studying the quantum correlations of 
Gaussian MPS of the form as in Eq. (|3}, and in relating 
them to the entanglement properties of the building block 7. 
The CM in Eq. Q describes a physical state if a; > 1 and 
s > Smin = (a; + l)/2 ll3ll . At fixed x, and so at fixed 
CM of mode 3 (output port), the entanglement in the CM 7^,^ 
of the first two modes (input port) is monotonically increas- 
ing as a function of s (as it can be checked by studying the 
respective symplectic eigenvalue rjss), ranging from the case 
s — Smin when 7^3 is separable to the limit s ^ 00 when the 
block 7^3 is infinitely entangled. Accordingly, the entangle- 
ment between each of the first two modes of 7 and the third 
one decreases with s. The main question we raise is how the 
initial entanglement in the building block 7 gets distributed 
in the Gaussian MPS F°"* . The answer will be that the more 
entanglement we prepare in the input port 7^^,, the longer the 
range of the quantum correlations in the output MPS wiU be. 
We start from the case of minimum s. 
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FIG. 2: (Color online) Entanglement distribution for a six-mode 
Gaussian MPS constructed from (a) infinitely entangled EPR bonds 
and (b) finitely entangled bonds given by two-mode squeezed states 
of the form Eq. J2j with r — 1.1. The entanglement thresholds 
with k = 1 (solid red line), k — 2 (dashed green line) and fc = 3 
(dotted blue line) are depicted as functions of the parameter x of the 
building block. For s > Sk, all pairs of sites i and j with j\ < k 
are entangled (see text for further details). 



Short-range correlations. — Let us consider a building block 
7 with s = Smin = (x + l)/2. It is straightforward to eval- 
uate, as a function of x, the Gaussian MPS in Eq. (|3} for an 
arbitrary number of modes (we omit the CM here, as no par- 
ticular insight can be drawn from the the explicit expressions 
of the covariances). By repeatedly applying the PPT criterion, 
one can analytically check that each reduced two-mode block 
7°"* is separable for |i — j| > 1, which means that the output 
MPS r°"* exhibits bipartite entanglement only between near- 
est neighbor modes, for any value of a; > 1 (for a; = 1 we 
obtain a product state). While this certainly entails that r°"* 
is genuinely multiparty entangled, due to the translational in- 
variance, it is interesting to observe that, without feeding en- 
tanglement in the input port 7^,^ of the original building block, 
the range of quantum coiTelations in the output MPS is mini- 
mum. The pairwise entanglement between nearest neighbors 
will naturally decrease with increasing number of modes, be- 
ing frustrated by the overall symmetry and by the intrinsic 
limitations on entanglement sharing (the so-called monogamy 
constraints We can study the asymptotic scaling of this 

entanglement in the limit x 00. One finds that the corre- 
sponding symplectic eigenvalue rji.i+i is equal to (N — 2) /N 
for even N, and [{N — 2) /N] for odd N: neighboring sites 
are thus considerably more entangled if the ring size is even- 
numbered. Such frustration effect on entanglement in odd- 
sized rings, akeady devised in a similar context in Ref. fl^, 
is quite puzzling. An explanation may follow from counting 
arguments applied to the number of parameters (which are re- 
lated to the degree of pairwise entanglement) characterizing a 
generic pure state on harmonic lattices |20]. 

Medium-range correlations. — The connection between input 
entanglement and output correlation length can be investi- 
gated in detail considering a general 7 with s > s,„in. The 
MPS CM in Eq. (|3} can still be worked out analytically for 
a low number of modes, and numerically for higher N. Let 
us keep the parameter x fixed; we find that with increasing s 
the coiTelations extend smoothly to distant modes. A series of 
thresholds can be found such that for s > Sk, two given 
modes i and j with |i — j| < fc are entangled. While trivially 
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FIG. 3; (Color online) Entanglement of formation between two sites 
i and j in a six-mode Gaussian MPS, with |i — i| equal to: (a) 1, (b) 
2, and (c) 3, as a function of the parameters x and d — s — Smin. 

si{x) = Smin for any N (notice that nearest neighbors are 
entangled also for s — si), the entanglement boundaries for 
k > 1 are in general different functions of x, depending on 
the number of modes. We observe however a certain regular- 
ity in the process: Sk{x, N) always increases with the integer 
k. These considerations follow from analytic calculations on 
up to ten-modes MPS, and we can infer them to hold true for 
higher N as well, given the overall scaling structure of the 
MPS construction process. Very remarkably, this means that 
the maximum range of bipartite entanglement between two 
modes, i.e. the maximum distribution of multipartite entangle- 
ment, in a Gaussian MPS on a translationally invariant ring, 
is monotonically related to the amount of entanglement in the 
reduced two-mode input port of the building block. 

To clearly demonstrate this intriguing connection, let us 
consider the example of a Gaussian MPS with iV = 6 modes. 
In a six-site translationally invariant ring, each mode can be 
correlated with another being at most 3 sites away (fc = 
1, 2, 3). From a generic building block Eq. Q, the 12 x 12 
CM Eq. (|3ji can be analytically computed as a function of 
s and x. We can construct the reduced CMs 7°"'i_j. of two 
modes with distance fc, and evaluate for each fc the respec- 
tive symplectic eigenvalue ?7i.i+fe of the corresponding partial 
transpose. The entanglement condition s > Sk will corre- 
spond to the inequality rii,i+k < 1- With this conditions one 
finds that S2(x) is the only acceptable solution to the equa- 
tion: 72s« - 12{x^ + 1)56 + {-Mx* + 28x^ - 34)3* + {x^ - 
5x* - 5a;2 + l)s'^ + (x^ - lf{x* - Qx"^ + 1) = 0, while 
for the next-next-nearest neighbors threshold one has simply 
s^{x) = X. This enables us to classify the entanglement distri- 
bution and, more specifically, to observe the interaction scale 
in the MPS r°"*: Fig. |2ja) clearly shows how, by increas- 
ing initial entanglement in 7^^, one can gradually switch on 
quantum correlations between more and more distant sites. 

We can also study entanglement quantitatively. Fig. |3l 
shows the entanglement of formation Ep of 7°"+^ for k — 
1,2,3 (being computable in such symmetric two-mode reduc- 
tions), as a function of x and d = s — Smin- For any (x, d) 
the entanglement is a decreasing function of the integer fc, 
i.e. quite naturally it is always stronger for closer sites. How- 
ever, in the limit of high d (or, equivalently, high s), the three 
surfaces become close to each other. We want now to deal 
exactly with this limit, for a generic number of modes. 
Long-range correlations. — In the limit s 00, the ex- 
pressions greatly simplify and we obtain a A^-mode Gaussian 
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MPS r°"* of the form Eq. gli, where C and C'^ are com- 
pletely degenerate circulant matrices, with (C^^)i i — — 
[{N - 1) + x^]/iNx), (C-i),,,^, = c, = -'l)/(iVx); 
and accordingly {C)ij = = [1 + (TV - l)a;2]/(iVa;), 
(C^^)i.j^i = Cp = —Cq . For any A^, thus, each individ- 
ual mode is equally entangled with any other, no matter how 
distant they are. 

The asymptotic limit of our analysis shows then that an in- 
finitely entangled input port of the building block results in a 
MPS with maximum pairwise entanglement length. These iV- 
mode Gaussian states are well-known as useful resources for 
multiparty CV communication protocols (l^l- The CM r°"* 
of these MPS can in fact be put, by local symplectic (unitary) 
operations, in a standard form parametrized by the single- 
mode purity nioc = [ciqap)^'^/'^ |9|]. Remarkably, in the limit 
IJ-ioc — *■ (i.e. X oo), the entropy of any K-sized {K < N) 
sub-block of the ring, quantifying entanglement between K 
modes and the remaining N — K, is infinite [9]. Within the 
MPS framework, we also understand the peculiar "promis- 
cuous" entanglement sharing |18] of these fully symmetric 
states: being them built by a symmetric distribution of infinite 
pairwise entanglement among multiple modes, they achieve 
maximum genuine multiparty entanglement while keeping the 
strongest possible bipartite one in any pair. Let us note that in 
the field limit (TV oo) each single pair of modes is in a 
separable state, as they have to mediate a genuine multipartite 
entanglement distributed among all the infinite modes fl. 

Finitely entangled bonds and experimental feasibility. — Vv^e 
finally discuss how to implement the recipe of Fig. ^ to pro- 
duce Gaussian MPS experimentally. The three-mode building 
blocks can be engineered for any choice of (x, s) (within the 
practical limitation of a finite available degree of squeezing) 
by combining three single-mode squeezed states through a se- 
quence of up to three beam splitters . The EPR measure- 
ments are in turn realized by homodyne detections fs*]. 

The only unfeasible part of the scheme is constituted by 
the ancillary EPR pairs. But are infinitely entangled bonds 
truly necessary? One could consider a F™ given by the di- 
rect sum of two-mode squeezed states of Eq. (|2j, but with fi- 
nite r. Repeating our analysis to investigate the entanglement 
properties of the resulting Gaussian MPS with finitely entan- 
gled bonds, we find that, at fixed (x, s), the entanglement in 
the various partitions is degraded as r decreases, as somehow 
expected. Crucially, this does not affect the connection be- 
tween input entanglement and output correlation length. Nu- 
merical investigations show that, while the thresholds Sk for 
the onset of entanglement between distant pairs are quantita- 
tively modified - a bigger s is required at a given x to com- 
pensate the less entangled bonds - the overall structure stays 
untouched. As an example, Fig.|2jb) depicts the entanglement 
distribution in six-mode MPS obtained from finitely entangled 



bonds with r = 1.1, corresponding to « 6.6 dB of squeezing. 
Single-bonded Gaussian MPS, which surprisingly encompass 
a broad class of physically relevant multimode states, can thus 
be experimentally produced, and their entanglement distribu- 
tion can be precisely engineered starting from the parameters 
of a simple bisymmetric three-mode building block, with the 
supply of Vwo-vaode, finitely squeezed states. 

Concluding remarks. — We have shown that the range of 
pairwise quantum correlations in translationally invariant iV- 
mode Gaussian MPS is determined by the entanglement in the 
input port of the building block. As a consequence of this 
interesting connection, a striking difference between finite- 
dimensional and infinite-dimensional MPS is unveiled, as the 
former are by construction slightly entangled for a low dimen- 
sionality of the bonds 1 10], and their entanglement is short- 
ranged 1 1 1 ] . We proved instead that pure, fully symmetric, iV- 
mode Gaussian states are exactly MPS with minimum bond 
cardinality {M — 1): yet, their entanglement can diverge 
across any global splitting of the modes, and their pairwise 
quantum correlations have maximum range. How this feature 
connects with the potential validity of an area law for critical 
bosonic systems is currently an open question. 
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